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A simple description for obtaining the parity distribution of nuclear levels in the pf 
+ shell as a function of excitation energy was recently derived. We implement 
this in a global nuclear level density model. In the framework of the statistical 
model, cross sections and astrophysical reaction rates are calculated in the Fe 
region and compared to rates obtained with the common assumption of an equal 
distribution of parities. We find considerable differences, especially for reactions 
involving particles in the exit channel. 



1. Introduction 

The nuclear level density is an important ingredient in the prediction of 
nuclear reaction rates in astrophysics. So far, all theoretical, global cal- 
culations of astrophysical rates assume an equal distribution of the state 
parities at all energies. It is obvious that this assumption is not valid at 
low excitation energies of a nucleus. However, a globally applicable recipe 
was lacking. We combine a formula for the energy-dependent parity distri- 
bution with a microscopic- macroscopic nuclear level density^ . The formula 
describes well the transition from low excitation energies where a single 
parity dominates to high excitations where the two densities are equal. It 
was tested against Monte Carlo shell model calculations. This treatment is 
further applied in the calculation of astrophysical reaction rates for nuclei 
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Figure 1. Left;. Ratio of odd to even parity partition function Z—jZj^ versus inverse 
temperature ,0 for ■''^Fc, ^"^Ni and "^Zn, calculated with tlie Monte Carlo Method^ .ijjg/ti: 
The parity ratio p_/p_)_ versus excitation energy E^. The solid lines are calculated using 
Eq. (1), and the solid circles are obtained with the Monte Carlo Method of NaJsada et. 
aP. 



in the Fe region. 



2. Method 

Y. Alhassid et. aP have calculated the ratio for nuclei in the iron 

region using the complete pf + /2 shell. In Fig. 1 the ratio of odd to even 
parity states as a function of inverse temperature /3 for three nuclei in the 
iron region is shown. The observed parity dependence can be explained 
quantitively by a simple model. Single particle levels are divided into two 
groups, according to their individual parities. The group which has the 
smaller average occupation number is denoted by 11. The distribution of 
the occupancies of the 11 parity group can be assumed to be Poisson, if the 
single particle states are occupied independently and randomly: 

P(n) = l-e-f . 
ni 

Here / is the average occupancy of orbitals with parity 11. The probability 
to have an odd/even parity state is given by: 

P+(n) = = e"^cosh/ 

n,even 

P_(n) = P{n) = e-^sinh/ , 

n,odd 
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and their ratio by 

^ = |^=tanh/. (1) 

The arguments leading to this equation are easily extended to the case 
where the protons and neutrons are treated separately — / has only to be 
replaced by the sum of individual contributions from neutrons and protons. 

The back-shifted Fermi-gas model (BBF) of Rauscher et. al^ is used 
to calculate the total partition function Z at constant excitation energy. 
Using + Z_ = Z and Eq. (1), wc can determine ^+/- and calculate 
the thermal energies for even- and odd parity states. Canonical entropies, 
heat capacities and the parity projected level densities are calculated from 
standard thermodynamic relations. 

3. Results and Discussion 

We have used the parity dependent level density to calculate astrophysical 
reaction rates in the global Hauser-Feshbach (HF) model NON-SMOKER^ 
and compared the results to the standard values^. We show the parity de- 
pendence for three nuclides, ^''Fe,^^Ni, ^®Zn and investigated all reaction 
channels involving these nuclei. The impact on the rates involving the Ni 
and Zn nuclei is small and negligible compared to the remaining uncertain- 
ties in the global HF model. This is due to the fact that a sufficiently large 
number of excited states is known experimentally. Up to 20 experimental 
states are considered in the standard calculation and only above the last 
known state, the theoretical level density is in effect. However, the case 
is (li£Fc;rc!nt for reactions involving *'^Fc. No information on experimental 
states is known here and therefore the full impact of the parity dependence 
can be seen. In the (n, 7) case 20% difference are found. Much larger differ- 
ences are seen in the reactions involving ^^Fe in the final particle channel. 
Because of lack of negative parities at low excitation energies, the particle 
emission channel becomes strongly enhanced in all such reactions with low 
or negative Q values. 

4. Outlook 

The case for ^"^Fe shows that a large effect of the parity dependence can be 
expected far from stability where no experimental information on excited 
states is available and that it is extremely important to include such a 
modified level density. The current approach is valid only for even-even 



February 8, 2008 11:14 WSPC/Trim Size: 9in x 6in for Proceedings 



cgsmocelj 



4 




Figure 2. Comparison of the parity dependent reaction rates to the standard®, which 
assumes an equal distribution of odd and even parity states. 



nuclei in the pf + gg /2 shell. Work is in progress to extend this description 
to be able to calculate the parity distribution for a large number of nuclei 
far from stability on the proton-rich as well as neutron-rich side. 
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